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We study the low energy edge states of bilayer graphene in a strong perpendicular magnetic field. 
Several possible simple boundaries geometries related to zigzag edges are considered. Tight-binding 
calculations reveal three types of edge state behaviors: weakly, strongly, and non-dispersive edge 
states. These three behaviors may all be understood within a continuum model, and related by non- 
linear transformations to the spectra of quantum Hall edge-states in a conventional two-dimensional 
electron system. In all cases, the edge states closest to zero energy include a hole-like edge state 
of one valley and a particle-like state of the other on the same edge, which may or may not cross 
depending on the boundary condition. Edge states with the same spin generically have anticrossings 
that complicate the spectra, but which may be understood within degenerate perturbation theory. 
The results demonstrate that the number of edge states crossing the Fermi level in clean, undoped 
bilayer graphene depends both on boundary conditions and the energies of the bulk states. 



I. INTRODUCTION AND PRINCIPAL 
RESULTS 

The integer quantized Hall effect is a generic behavior 
of two-dimensional electron systems in a strong perpen- 
dicular magnetic fielciS^. The primary manifestation of 
the effect is a precise quantization of the Hall conduc- 
tance axy to integer multiples of e^//i, with coefficient 
determined by the electron density. That this system 
can carry current at all is in some ways surprising, be- 
cause the spectrum of the bulk system takes the form of 
Landau levels, highly degenerate states at discrete energy 
values, with gaps separating these isolated sets of states. 
With chemical potential placed in any of these gaps one 
naively expects the system to be insulating. A basic 
explanation for the existence of Hall currents in clean 
systems involves edge stated the energies of the Lan- 
dau level states disperse as their guiding center quantum 
numbers- X approach the physical edge of the sample, 
and are thus current-carrying in a particular direction for 
a given edge. A difference in the occupation of states at 
opposite edges of the sample leads to a net current, with 
a voltage difference perpendicular to that curreniP, such 
that their ratio yields a Hall conductivity quantized at 
the number of distinct edge state branches which cross 
the Fermi level at a given edg^. When the chemical po- 
tential lies in an energy gap in the bulk of the system, 
the only low energy excitations of the system are present 
at its edges. These dominate the low-energy physics of 
the system. 

More recently, it has been recognized that the presence 
of gapless edge states in a system with a bulk energy 
gap is the defining characteristic of a more gen eral class 
of systems, known as topological insulatord^. Interest- 
ingly, in such systems states with different quantum num- 
bers at the same edge may cross the Fermi energy such 
that they carry current in opposing directions, so that 
there are both hole-like and particle-like currents at the 



same edge. These states can be topologically protected 
from backscattering, and allow the transport of currents 
without dissipation. In addition, in such systems one 
may observe transport of quantities other than electric 
charge (e.g., spin) along their edges while carrying no 
electric current. The realization of such currents would 
be major step in the exploitation of d egrees of freedom 
beyond charge in electronic device^^. 

One system known to possess this sort of behavior is 
graphene. Graphene is a two-dimensional honeycomb lat- 
tice of carbon atoms, which recently has become available 
in the laboratorjii^lMi, Electronic states near the Fermi 
energy in this system largely reside in orbitals of the 
carbon atoms, and when undoped, the low energy con- 
tinuum description of the electron states is best given in 
terms of the Dirac equation^'' With an appropriate 
spin-orbit coupling term, it was shown that single layer 
graphene could become a topological insulator even in 
the absence of a magnetic field. 1^. However, subsequent 
estimates of the strength of this spin-orbit coupling in 
real graphene suggested that the effect would be very 
difficult to observtPHUl. 

Crossing of edge states with different quantum num- 
bers can nevertheless be re alized in single layer graphene 
in the quantum Hall regim^i2'20| 'pj-^jg jg (j^g j^g unique 
Landau level spectrum, which has both positive and neg- 
ative energy states (and is particle- hole symmetric) , with 
the former supporting upwardly dispersing edge states 
and the latter downward dispersing edge states. In the 
absence of interactions and Zeeman coupling, there are 
four Landau levels precisely at zero energy in the bulk, 
with each spin state supporting a particle- like and a hole- 
like edge state at each edge. When Zeeman coupling is 
included, the two spin states split so that one hole-like 
state crosses one electron-like state at each edge. This 
allows for dissipationless spin transport at the edge^. 
The inclusion of electron-electron interactions transforms 
the crossing edge states into a magnetic domain wall with 
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FIG. 1: {Color online) Energy bands of a Bernal-stacked 
graphene bilayer nanoribbon, A'^ = 240 atoms across in each 
layer, each with zigzag edges. Perpendicular field is lOOT, 
interlayer bias V — 0.00ft, and 71 = 0.25t, with t the in- 
plane hopping amplitude. 



Luttinger liquid propertied^. Tliis structure may explain 
the presence of apparently metallic behavior for undoped 
graphene in magnetic fields of order ^ip^jUMSI l^ whi ch 
gives way to an insulating state in stronger fieldd^HHl, 

The rich physics associated with crossing edge states 
suggests that one may expect to find unusual behaviors 
in other systems that support both particle- and hole- 
like edge states. In this context, bilayer graphene is a 
particularly interesting candidate to investigate. Even 
in the presence of interlayer coupling, bilayer graphene 
supports (in the absence of Zeeman coupling) eight zero 
energy stated^!!. Unlike the single layer case, this degen- 
eracy can be broken and controlled using an external, 
perpendicular electric fielcP^. This raises the possibility 
of controlling the edge state structure via a combina- 
tion of this electric field and the Zeeman coupling (which 
may be manipulated using a parallel magnetic field). In 
what follows, we investigate the edge state structure of 
a bilayer graphene ribbon using both tight-binding cal- 
culations and the Dirac equation, assuming appropriate 
boundary conditions for the latter. We focus on ribbons 
with zigzag edges, as well as some simple extensions of 
this involving "bearded" edged^ll at a given edge. 

Fig. [T] illustrates a typical spectrum for the energy 
states of a graphene bilayer nanoribbon in a perpendic- 
ular magnetic field, as a function of ky, the wavevector 
along the ribbon. One may see two regions supporting 
very flat bands in the vicinity of the y component of the 
vectors K and K', the locations of the Dirac points in 
the Brillouin zone of a graphene sheet. These resul ts are 
consistent with those obtained in previous studied^SlMI 
The bulk states for the valley on the left appear at ener- 



gies El = V/2, where V is the potential energy difference 
between the layers due to a voltage bias, and, in the limit 
of small V, 
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Here 71 is the hopping amplitude between overlaid sites 
of the Bernal-stacked layers, and uJc = \/2Tivf in which 
I = \Jch/eB is the magnetic length associated with the 
perpendicular magnetic field B, and is the speed of 
electrons in the vicinity of a Dirac point in the absence 
of interlayer hopping. Analogous bulk energy states are 
present around Ky at e'l — — ei and e'2 = — £2- 

From the form of ei it is clear that the wavefunctions 
corresponding to this band reside in a single sheet. For 
the right edge one finds no dispersion in this energy band, 
so that this edge state cannot contribute to the Hall con- 
ductivity of the system. Such non-dispersive edge states 
are one type of behavior that is supported by the bilayer 
graphene edge, and are very analogous to those of the ze- 
roth Landau level of a single graphene layer with a zigzag 
edg^. 

At the same edge, £2 disperses downward, and we shall 
see that its dispersion has the approximate form 
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where /(AT), which may be determined variationally, 
grows monotonically from 1 when X is deep in the sys- 
tem bulk to large positive values when X is well over 
the system edge. This means that one expects the edge 
state to disperse downward, from the bulk value £2 to 
a value close to —V/2, as is apparent in Fig. [ij Be- 
cause the range of energies available to these edge states 
is limited, they disperse relatively slowly, and represent a 
second type of edge state that is supported by the bilayer 
graphene system. 

On the left edge of the system, there is an edge state 
which originates in the K' valley at —£2, and approaches 
V/2 as X moves well outside the bulk [in analogy with 
Eq. ([2])]. Rather than becoming degenerate with £1, this 
begins to disperse downward as the edge is approached, 
so that ultimately there are both particle-like and hole- 
like states dispersing from the vicinity of £1. This is 
analogous to the single layer casd^, for which a zigzag 
edge supports both particle-like and a hole-like branches 
dispersing from the n = Q Landau level. Note that these 
states disperse rapidly toward ±00 as the wavefunction 
centers move across the edge, representing a third type 
of behavior supported by this system, and is most simi- 
lar to behaviors apparent in conventional quantum Hall 
systems^'. We will see below that these states are most 
simply understood in terms of the 71 = single layer edge 
states, coupled together by 71, resulting in level repulsion 
and anticrossings. 

This complicated structure suggests interesting possi- 
bilities for the low-energy edge states in bilayer graphene. 
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For Fermi level precisely at zero energy {i/ = 0) and V 
exceeding the Zeeman splitting, one finds counterprop- 
agating edge states for each spin, one from each valley, 
at a given edge. This contrasts strongly with the ly ~ 
state of a conventional two-dimensional electron system, 
for which there are no edge states at all. In principle the 
counterpropagating states will mix and localize due to 
disorder, but because they are well-separated in ky, the 
localization length could be relatively long. Thus charge 
transport due to these edge states might be observable 
over short distance scales. It is also possib le that they 
could be observed in thermal transpori!22lMI 

On the other hand, for large Zeeman splitting Ez and 
small V, both electron-like states above {V — Ez)/2 for 
spin up states will cross the hole-like states below {Ez — 
V) /2 at zero energy. In the absence of perturbations 
that can admix different spin state^^H^ these channels will 
remain open so that v ^ would become a quantized spin 
Hall statei'^. 

Finally, it is interesting to note that if the ratio "fi/ujc 
can be tuned below 1, 62 would fall below 0, and no edge 
states would cross the Fermi level at all when v = ii 
the Zeeman coupling is sufhciently small. In principle 
this can be accomplished with large magnetic fields, but 
would require values well above those currently available 
in the laboratory for the bare value of 71. It is possible, 
however, that the effective value of 71 could be decreased 
by an in-plane magnetic field. Presuming the energy 62 
can be made to cross through zero energy for the un- 
doped system, this leads to the possibility of driving a 
topological phase transition within the v = Q state. The 
change in the edge state structure for such a transition 
would be accompanied a bulk change in the state, from 
a partially valley-polarized to an unpolarized state. 

The remainder of this article is organized as follows. 
In Section [n] we describe our tight-binding results for the 
edge state structure in more detail, and show how the 
results evolve from the single layer results'^^ as interlayer 
tunneling is turned on from zero. Section [lIT| discusses the 
continuum representation of these results. We conclude 
with a summary and some speculations in Section [Tvl 



II. NUMERICAL RESULTS FOR THE 
TIGHT-BINDING MODEL 




FIG. 2: (Color online) Unit cell for Bernal-stacked graphene 
bilayer nanoribbon, with zigzag edges. Dashed lines indicate 
bonds on lower layer, solid lines are bonds on upper layer. 



74, whose effect we assume to be negligible in the presence 
of a perpendicular magnetic fielcpxl. We consider the unit 
cell structure, which has width a, to be infinitely repeated 
in the y direction, and to be repeated a finite number 
of times in the x direction. The resulting structure has 
zigzag edges in both layers on both sides of the ribbon. 
Other edge constructions can be generated by removing 
atoms at the edge from the top or bottom layer. Remov- 
ing an odd number of atoms from one of the layers in 
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this way generates a "bearded" edge^ removing an even 
number returns the edge to a zigzag form. We explore 
two such constructions below. To implement the mag- 
netic field, we introduce a vector potential into the hop- 
ping matrix element between neighboring atoms a and h 
in the standard way, t — )• t exp A • dr], where A is 

the vector potential associated with the magnetic field, 
and we have taken h = \. Note that in order to avoid 
using excessively large numbers of atoms in a unit cell, 
we set the magnetic field to be rather large {B =100T), 
so that our ribbon is several magnetic lengths across. Al- 
though this is beyond what is typically attainable in the 
lab, our results should be qualitatively the same as for 
wider ribbons in lower magnetic fields. 



A. Zigzag Edges 



Our numerical calculations are based on a simple near- 
est neighbor tight-binding model for graphene, with hop- 
ping amplitude t which we take as our unit of energy in 
what follows. The basic unit for the bilayer crystal struc- 
ture is illustrated in Fig. [2] in which there is an upper 
and lower layer whose bonding structure is indicated. In 
addition there is a hopping matrix element 71 connecting 
sites lying above/below one another [red {A) and yellow 
{B) atoms in Fig. [2]. The graphene bilayer may also 
have longer range interlayer hopping parameters 73 and 
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FIG. 3: {Color online) Detail of tight-binding energy bands 
of zigzag bilayer ribbon, near left edge of system, for state 
emerging from Landau levels of the K valley. Unit cell con- 
tains 480 atoms, perpendicular field is lOOT, and V = 0.001 
(in units of t.) Results illustrated for several values of 71. 



Given the form of the tight-binding model, it is clear 
that there should be a continuous evolution of the spectra 
from that of decoupled layers (71 — 0) to the form exhib- 
ited in Fig. [1] for physical values of 71 . Fig. [3] shows an 
example of this for a series of 71 values, from 71 = 0.05 1 
to 71 = 0.2 1. Note that guiding center coordinates for 
the single particle states connected to the bulk states at 
K have the form X = {ky — Ky)£'^, up to an overall con- 
stant. Where the bands begin to strongly diverge from 
their bulk energies as a function of ky, the guiding cen- 
ter coordinate comes close to the physical edge of the 
system. This is easily be confirmed by the form of the 
wavef unctions. 



For the smallest value [Fig. [Sj^a)], it is clear that the 
basic structure of the spectrum involves particle-like and 
hole-like edge states, each dispersing from bulk bands 
around e — ±V/2. The two states converging towards 
zero energy are admixed by 71, creating an anticrossing. 
Note the gap associated with this anticrossing is rela- 
tively large, because 71 >> V. Thus one sees the spec- 
trum is largely similar to that of two uncoupled layers at 
different constant potentials. For all the results shown 
in Fig. [3] when X is sufficiently inside the bulk that the 
effect of the edge is quite small, one may see that the 
two levels closest to zero always initially approach one 
another as X moves towards the edge. The two modes 
then anticross, and furthermore anticross with the levels 
closest to ±.V/2. Interestingly, the two modes at ±.V/2 
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persist to slightly larger values of ky before diverging to 
large values of |e|. Note that of these two modes, the 
positive energy one is an edge mode of the bulk band in 
the K valley at e — V/2, while the negative one is the 
continuation of an edge state associated with a bulk band 
at - V/2 for the K' vahey. 

Edge states associated with the K valley for the other 
side of the ribbon behave relatively smoothly compared 
to the above, and are plainly visible in Fig. [T] This con- 
sists of a dispersionless edge state at e = V/2 associated 
with the bulk band £i, and an edge state dispersing down- 
ward from £2 toward —V/2, where it continuously joins 
to the particle-like branch of the edge states (for 71 — >• 0) 
associated with the bulk state at —V/2 of the K' valley. 
As we discuss in Sec. |IIH the behavior of these two states 
can be understood in a relatively straightforward man- 
ner from the continuum description of this system with 
appropriate boundary conditions. 



B. Variants on the Bilayer Zigzag Edge 

We next discuss the edge state spectra for two variants 
of the zigzag edge, known as the "bearded" edg^^Sl, This 
structure is created from a zigzag ribbon edge by remov- 
ing the outermost atoms at the edge. The structure can 
also be created by adding single atoms to the outermost 
points of zigzag edge. 

In the bilayer, structures involving bearded edges nat- 
urally emerge if one cuts all the bonds along a line in 
the zigzag direction. In addition to the zigzag geometry 
illustrated in Fig. [T] two other possibilities arise, as illus- 
trated in the insets of Figs. |4]and[5] The edges in these 
two latter cases both involve a single zigzag edge in one 
layer, and a bearded edge in the other. Unlike the rib- 
bon with two zigzag edges, these ribbons present atoms 
on the same sublattice at both edges. The difference be- 
tween the two zigzag-bearded edge ribbons is that in one 
case the atoms at an edge are uncoupled between layers, 
whereas in the other case the two outermost atoms form 
an interlayer dimer. 

The spectrum of the former case is illustrated in the 
main panel of Fig. [4] Prominently visible are bands of 
constant energy precisely at ztV/2. Such bands across 
the Brillouin zone are also visible when 71 = 0, the spec- 
tra of two single layer ribbons at potentials ztV/2 each 
with one bearded edge and one standard zigzag edge (see 
Fig. [6]). In terms of a continuum model, this latter re- 
sult has a simple interpretation: for the K' valley of the 
bottom (—V/2) layer, the boundary condition may be 
taken to be vanishing of the A sublattice component on 
both edges, leading to dispersionless edge states on both 
sides. In this structure the dispersionless state of the 
left edge continues through the K valley, where it has no 
simple continuum interpretation in terms K valley states. 
These states are very localized on the edge atoms of the 
bearded edge, and because their hybridization with the 
rest of the ribbon is extremely weak, and there is no 



hopping directly among them, the energy of the state is 
essentially pinned at —V/2. 

For the K valley, the boundary condition in the same 
layer is 5=0, so that one finds the pair of dispersing 
particle-like and hole-like edge states of the zeroth Lan- 
dau level for a standard zigzag edg^^at both edges. Note 
the unusual situation that three bands are degenerate at 
—V/2 near the Ky point; the extra state is most natu- 
rally interpreted as a continuation of the n — Landau 
level edge state from the K' valley. 

This situation evolves in a simple way when 71 is in- 
creased from zero. For the K valley, the bulk mode at 
£1 = V/2 is localized on a single sublattice which is not 
directly affected by the boundary conditions, and so re- 
mains dispersionless at both edges. The other two K 
valley levels which were degenerate at —V/2 for 71 = 
evolve into a bulk mode at 82, which has particle- like 
edge states due to the boundary condition, and into an 
edge mode whose energy remains near —V/2 for ky suf- 
ficiently close to Ky, but develops a strong hole-like dis- 
persion away from the valley center. One also observes 
the edge state from the K' valley at —V/2. 

It is interesting to contrast this edge state structure 
with what is apparent in Fig. [T] In addition to being 
considerably simpler, the edge state structure of Fig. [4] 
has no slowly-dispersing edge states, as is the case for the 
other edge constructions we consider. Moreover, there 
are no edge states of any kind crossing the Fermi level 
when it is at zero energy in this particular case. This 
demonstrates that in bilayer graphene, one may or may 
not have edge states crossing zero energy for the same 
bulk spectrum, depending on boundary conditions. In 
the former case these are counterpropagating, so that no 
charge current is present at the edge in equilibrium, al- 
though these may transport energjEH, That the presence 
or absence of low-energy edge excitations can depend on 
boundary conditions is somewhat unusual for a quan- 
tum Hall state, but is allowed because there are no strict 
quantum numbers distinguishing the counterpropagating 
states. When counterpropagating edge states carry dif- 
ferent quantum numbe rs (e.g . , spin) we expect their pres- 
ence to be more robustPiESl 

Finally, we consider the situation in which the outer- 
most atoms at the edge are dimers, tunnel-coupled by 71. 
The corresponding spectrum is illustrated in Fig. [5] In 
this situation there are no dispersionless states because 
the boundary conditions involve the sublattices on which 
the bulk states at £1,2 for the K valley (and —£1,2 for the 
K' valley) reside. Interestingly, we find two edge states 
that "thread" the gaps between the bulk states. Unlike 
the previous case, where each extra atom of a beard con- 
nected to atoms only through a single bond, in this case 
these atoms are coupled to the zigzag edge of the op- 
posing layer through 71. Thus it is not surprising that 
states localized on these sites would develop a dispersion, 
whereas in the previous case there was none. This situa- 
tion is rather unique in supporting quasi-one dimensional 
states at the edge which are not directly connected to any 
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FIG. 4: {Color online) Energy states near zero for graphene 
ribbon with one layer bearded at each edge. Unit cell contains 
474 atoms. Perpendicular field is lOOT, V = 0.001, 71 = 0.25 
(in units oft). Upper inset illustrates edges of unit cell. 
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FIG. 6: {Color online) Energy states near zero for graphene 
ribbon with one layer bearded at each edge. Parameters are 
identical to those of Fig. |4] except 71 = 0. 



III. CONTINUUM DESCRIPTION 
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FIG. 5: {Color online) Energy states near zero for graphene 
ribbon with one layer bearded at each edge. Dimer atoms 
protrude at edges in this construction. Unit cell contains 478 
atoms. Perpendicular field is lOOT, V = 0.001, 71 = 0.25 (in 
units of t). Upper inset illustrates edges of unit cell. 



bulk state. 

The behavior of the dispersive energy levels in each 
of the above mentioned edge structures (Figs, [ij [4] and 
[5]) can be understood within a continuum theory with 
the appropriate boundary condition. This is described in 
detail in the next section. 



We consider a Bernal-stacked bilayer graphene ribbon 
of finite width L in the i-direction, where inter-layer hop- 
ping is assumed to be only between the overlaid sites [red 
{A) and yellow {B) in Fig. ^ with an amplitude 71 , and 
an inter-layer voltage bias V is appHed. Using a basis of 
4-spinors {B,A,B,A), where A, B {A, B) denote wave- 
function components on sublattices A and B of the top 
(bottom) layer, the Dirac Hamiltonian projected onto a 
given ky in the vicinity of the valley K is given by the 
4 X 4-matrix 



/ -V/2 \ 

LJco) -V/2 71 

71 V/2 u^a 

V u^a^ V/2) 



(3) 



Here a ~ 



(x-X)] and at = + (.t - 

X)], where x and the guiding-center coordinate X = 
i{ky — Ky) are in units of the magnetic length i, and 
(jJc — \f2'hvp/l. In the vicinity of the other valley (K' 
point), the same Hamiltonian (with V — > —V) applies in 
the basis of inverted 4-spinors (A, B, A, B). As already 
discussed in Sec. |lj for ^ ^ TijWc the bulk solution 
for the energy spectrum of ([S]) includes two low energy 
levels. El = V/2 and £2 [ Eg. ([T])]. The corresponding 
eigenf unctions are given bjESl 



l*i> = 




1*5 
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in which \n) = <j)n{x — X) are the harmonic oscillator 
wave-functions, and A/" is a normalization factor. The 
dispersion of when X approaches the edge can be 
found by imposing the appropriate boundary condition 
aX X — ±L/2. Below we study separately four distinct 
boundary conditions, compatible with the tight-binding 
calculations of previous section. 



A. Right Zigzag Edge: B{L/2) = B{L/2) = 

In a bilayer ribbon with zigzag edges including an 
integer multiple of unit cells, the boundary conditions 
at the right and left edges are fundamentally different. 
We first consider the right-hand edge [x — i/2), at 
which the wave-function is forced to vanish on the B 
sublattice of both layers. We therefore look for solu- 
tions of the form '^f3''{x) = {B{x),A{x),B{x),A{x)) 
where B{L/2) = B{L/2) = 0. From Eq. Q it is obvi- 
ous that the bulk wave- function 'S'l already obeys this 
boundary condition, hence £i is non-dispersive in anal- 
ogy with the monolayer case. In contrast, the component 
B{x) of ^2 is non- vanishing; however it is smaller than 
the components A{x), A{x) in the small V limit. This 
suggests that ^^'^^'^ is given by a smooth deformation 

of '4'2, which dictates a dispersion 62'^^'^ {X) of the cor- 
responding eigenvalue. An exact analytic evaluation of 
^edge^^-j is uot possible. However, as we show next, an 
approximation based on either a variational calculation 
or a perturbation expansion in the inter-layer hopping 
can explain the right-hand dispersion of £2 in Fig. [l] 

We start with a variational approach, similar to the 
one adapted in Ref. [3T]for a single layer graphene. The 
variational ansatz on ^2'^^'^{x) is taken to be the simplest 
modification of the bulk function '$'2(2:) which obeys the 
boundary condition. We therefore assume B{x) = 0, and 
apply the variational principle to the remaining three 
components, out of which only B(x) is restricted by 
the vanishing boundary condition. Note that since the 
spectrum of the Dirac Hamiltonian is unbounded, the 
standard procedure of minimizing the energy expecta- 
tion value £ = (H) is not applicable. However, it turns 
out possible to express it as a monotonic function of an 
"effective energy" functional with a well-defined mini- 
mum. To see this, we first impose the extremum condi- 
tion Se/6A = Se/SA = which yield 



where 



A{x) 



71 



{e + V/2) 



B{x) , 



(5) 
(6) 



Evaluating (H) for this state, in the small V limit, pro- 
duces an expression for £ as a functional of B{x) only: 



y (7?-c.g(««^)B) 
2 (7?+c.2(aat)s) 



(7) 



{aa^)i 



J dxB*{x)aa^B{x) 
JdxlBixW 



(8) 



Quite interestingly, the expectation value {a^a)B (implic- 
itly dependent on X via the definition of a, a^) is equiv- 
alent (up to an additive constant) to the energy of a 
quantum Hall edge states in an ordinary 2D electron gas. 
In particular, it is identical to the functional associated 
with the square of the energy of edge states in single layer 
graphene"^ ^, and can be minimized using a standard vari- 
ational ansatz for B{x). Notice that minimizing (a^a)^ 
with respect to B also minimizes £^ in Eq. [t] giving esti- 
mates for the states closest to zero energy. The dispersion 
curve f{X) = 1-1- min{{a^ a) b} has a known qualitative 
behavior as a function of X: in the bulk, B{x) = |0) 
hence f{X) = 1; as X approaches the boundary, f{X) 
increases monotonically and acquires large positive val- 
ues when X is well beyond the edge. When substituted 
in Eq. (|8|, this yields the dispersive energy band 



^edge 



V {^l-u:lJ{X)) 
2 {^l+u^lf{X)) 



(9) 



which decreases monotonically with X from the bulk 
value £2 to the saturated value e'^^'^{X) — —V/2 as 
/(X)^oo. 

An alternative approach to the derivation of the above 
dispersion law involves a perturbative expansion in 71. 
This approach turns out useful to develop insight about 
the prominent qualitative features of the spectrum for 
more complicated boundary conditions as well, even in 
the regime where it is not strictly justified to assume 71 
small. To this end, we define 



0\ 

71 
71 

V 00/ 



(10) 



as a perturbation on Hg = ^71=0 describing the uncou- 
pled layers. The eigenstates of Hq are single-layer Lan- 
dau level (LL) states. Focusing first on bulk states, the 
zero LL states and corresponding energies (split by the 
inter-layer bias V) are given by 



I* 



(0)\ 




-(0) 



I* 



(0)\ 




■^2 



V 

'2 



V 



(11) 



(12) 
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Since Iit\"i'f'^) = 0, the perturbation does not couple the 



top layer state [Eq. ( 11 )] to higher LL's, so that l^'i"'') = 
l^-i) [Eq. ^] and remains fixed at ei = V/2 for 
arbitrarily large 71. In contrast, if couples the bottom 
layer state [Eq. (12|] to the n — ±1 LL states in the top 
layer 



l^±i) 




V . 

e±i = ± 



(13) 



To second order in perturbation theory (and leading or- 
der in V/ujc), the resulting correction to S2 is 



-(2) 



E 

n=±l 



|(<I>„|ffi|$o)|^ 7?^ 



-(0) 



the resulting 



-(0) 

'-2 



(14) 



(2) 



To leading order in ji/uj, 
coincides with Eq. Q. 

We next consider edge states where X approaches the 
right edge boundary L/2. Since both |$o) and |(f>o) [Eqs. 
(11), (12|] have vanishing components on the B sublat- 



tice, the boundary condition is obeyed and e^i \ Sj^"* do 
not disperse. However, higher LL states are modified and 
consequently so is the energy eigenvalue £2 at finite 71. 
For 7i = 0, the wave-functions and energies ( 13 1 become 



perturbation expansion parameter becomes increasingly 
smaller, i.e., the coupling between layers effectively weak- 
ens. This behavior turns out to be valid for all types of 
boundary conditions. In the present case, we conclude 
that the dispersion curve S2{X) is monotonically decreas- 
ing and asymptotically approaches —V/2 for X — > 00, in 
agreement with the variational result Eq. ([9]). 



B. Left Zigzag Edge: A{-L/2) = A{-L/2) = 

The boundary condition on the left edge of the rib- 
bon, A{—L/2) = A{—L/2) = 0, creates a much stronger 
disturbance for both electronic wavefunctions [^'2) 
when X is close or to the left of —L/2, and changes 
their shape significantly. To analyze this case, we imple- 
ment the perturbative approach introduced in the previ- 
ous subsection. The uncoupled layers states |<&o), j^o) 



and the corresponding energy levels [Eqs. (11), (12)] are 
now split into two branches each: 



1$ 



±0/ 



1 



( 



\ 



± 



v 



|0l 



) ) 



V 



o(X)--±a.,yA(X), 



(18) 



1 

Mr 



-±i 



/ 



^±v/l + A(X)o 



(15) 



where \Qr)\x=l/2 = so that a^a\QB) — \{X)\Qb) with 
\{X) > the dispersion curve of a conventional lowest 
LL edge state, |lij) = o^IOr), and Mr is a normaliza- 
tion factor. Neglecting the contribution of higher LL, we 
obtain the second order correction to £2 



-fl(2) 



a.2[l + A(X)] 



(16) 



Note that Eq. ( 16 ) is similar to ( 14 ), with the expansion 



parameter 71 /ujc replaced by the X-dependent parameter 
7i/(Iic(X), where 



OJc{X) 



Nruj,^1 + \{X) 
^l(Oi?|0)| 



(17) 



When X is pushed farther towards the edge, u)c{X) is 
monotonically increasing due to a combination of the in- 
crease of \{X) in the numerator and the suppression of 
the overlap |(Ofl|0)| in the denominator. For X far be- 
yond the physical edge, Cjc{X) — 00. Hence, even in the 
physically relevant case where ^i/lOq > 1, the effective 



1 

Mo 



\ 



|0l; 





^±0 



V 



{X)^--±u:,./X{X) 



(19) 



where |0L)|,^_i/2 = 0, aU\{)L) - A(X)|Ol) [with \{X) 
the same as \{—X) of Eq. (15)] and = <A^l) a wave- 
function strongly confined to the edge. Note that the 
hole-like dispersive branch of the top layer state [£^q(X)] 
and the particle-like branch of the bottom layer [£:^q(X)] 
cross at zero energy. When we next turn on a finite 
but small inter- layer hopping 71 , these two branches mix 
and a gap will open up, yielding an avoided crossing as 
observed in Fig. [sja). For larger 71, each of the mix- 
ing branches separately will get modified and the band 
structure becomes more complicated. To leading order 
in perturbation theory, we consider the corrections due 
to mixing with higher LL states 



1*^1) 



1 



± 



V 



.±i--±Wc\/Ai(X) 



(20) 
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l^ii) 



1 

4i 



|1'l> 




V 



(21) 



where |l'i)U=_L/2 = 0, o)a\V^) = Ai(X)|l'i) with 
Ai(X) > 1, and = ajl'^). This yields the fol- 

lowing approximations for the hole-like and particle-like 
branches dispersing from the bulk energy levels ei, 62- 



V , 27f|(l'J4)p(±^cyA(X) + T/ 

V , 27?|(0'jOi)p(±..,yA(X)-l/ 

± i^c^\{X) ^ ^ 

2 NlMl\^{X) (c.2(Ai(X) - \{X)) ± 21/^,yA(X) 



(22) 



(23) 



In particular, the hole-like branch £\ _(X) and the 
particle-like branch e\^{X) develop a non-trivial (possi- 
bly non-monotonic) dependence on X, which shift their 
crossing away from zero energy. The gap opening at the 
avoided crossing point is given to leading order by degen- 
erate perturbation theory as 



Al(^) 



71 

v/A(X)AA2 



I(OlKl)| 



(24) 



As 71 becomes bigger, the second order corrections in 
Eqs. (22 1, (23 1 become increasingly dominant, and in 



particular the negative correction to £2 +(^) '^^^ \e^sA 
to the features observable in the spectrum depicted in 
Fig. [3]jb)-(d). However, it should be noted that (as in 
the previous case of boundary conditions, and for the 
same reason) the perturbative expansion systematically 
improves for the farthest edge states (corresponding to 
X very close to or beyond the left edge). The lowest 
energy levels are then approxim ated by the particle-hole 
symmetric values ±(y/2 — Wc\/A(X)), consistent with 
Fig. [31 



C. Top-Layer Bearded Edges 

B(±L/2) = i(±L/2) = 

The next type of boundary condition corresponds to 
the edges depicted in Fig. [4j A special feature of this 
particular configuration is that an identical (vanishing) 
boundary condition is imposed on both wave-function 
components associated with the overlaid sites of the inter- 
layer dimer, i.e. B{±L/2) ^ A{±L/2) ^ 0. Therefore, 
one can find a consistent solution to the Dirac equation 



= el*) [where H is given by Eq. ([s])] with A{x), 
B{x) being given by the same function (up to a constant 
prefactor). 

To see this, we note that the Dirac equation can be 
cast as a set of four coupled equations: 



uJco} B + 71 _B 
jiA + LOcaA 

UJrO^ B 



(e + y/2) B 
{e + y/2) A 
(e - F/2) B 
(e - V/2) A 



(25) 
(26) 
(27) 
(28) 



which can be combined to yield two coupled Schrodinger 
equations for the components A, B: 

{ujla)a-{e + V/2f)A = --ii{e + V/2)B (29) 
{uilaa^ -{e-V/2f)B = --^i{e - V /2)A . (30) 



Clearly, there is a solution to Eqs. (29), (|30j) of the form 
A — c^|Oe), B = csjOe) in which ca, cb are constants 
and |0e) is an eigenstate of the operator a satisfying the 
boundary condition. For X close to (or beyond) one of 
the edges ±L/2, the function |0e) satisfies the boundary 
condition |0e)|x=±L/2 = and the Schrodinger equation 



a^a|Oe) = A(X)|Oe 



(31) 



here A(X) is the same dispersion curve introduced in the 
previous subsections, corresponding to the edge disper- 
sion of a conventional LLL edge state. In fact, jOg) co- 
incides with jOij) (section [III A ) for AT > 0, and |0l) 
(section 
Eqs 



niB| 



for A < 0. Substituting this ansatz in 
(|29|) and ((30|, we get an eigenvalue equation for e: 



{e-V/2f~u:l[l + \{X)] 
li{e-V/2) 



li{s + V/2) 



{e + V/2Y-ulX{X) 



(32) 
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For V <C i^c,7i, the two lowest energy solutions are 



1 



r2 + 1 + 2X{x) 



-|±Wy2 fll_xix)[i + xix)] 



+ c^2A(X)[l + X{X)] [r2 + 1 + 2X{X)] 



(33) 



where F = 71 /wc- 

We first note that the above calculation recovers the 
known bulk solution for X{X) ~ and |0e) — |0). Indeed, 
Eq. ([33]) then yields e+{X) = £2 [Eq. The appar- 

ent second solution £_ ~ —V/2 does not correspond to 
a valid solution of the original Di ra c eq uation: inserting 
A = ca\0), B = cb|0) in Eqs. (25), (26) gives an ambigu- 
ous expression for the B component. (This can be traced 
back to an assumption that a)a\Qe) 7^ 0, which is not the 
case when |0e) is a bulk lowest Landau level state.) We 
therefore conclude that s±{X) converge to a single bulk 
energy level £2, which (as noted earlier) has evolved from 
the zero Landau level bulk state of the uncoupled bottom 
layer. However, as soon as X{X) is finite, Eq. (33) dic- 



tates that the bulk state splits into two dispersive bands: 
£_i_(X) is particle-like, and steeply deviates upward from 
£2 as X approaches the edge, i.e. with increasing X{X); 
£-{X) is hole-like, and steeply deviates downward from 
—V/2 as X{X) increases. This behavior is clearly seen in 
Fig. m 

We finally comment that in addition to the above men- 
tioned dispersive energy bands, there exists a trivial so- 
lution to this boundary problem where A ^ B ~ Q. 
Similarly to the case discussed in section III A this cor- 
responds to the bulk wave-function [see Eq. ([4])] 
which is not affected by the boundary. As a conse- 
quence, there is no dispersion of the bulk energy level 
£1 and it is maintained fixed at V/2 for arbitrarily large 
The other valley (K' point) contributes another 
non-dispersive state at energy —V/2, which corresponds 
to an eigenfunction localized on the B component 
only. Together with e±{X), this explains the entire spec- 
trum depicted in Fig. |4] 



D. Bottom-Layer Bearded Edges 

A{±L/2) = B{±L/2) = 

The boundary condition corresponding to the edge 
structure depicted in Fig. [sjcan be cast as A{±L/2) = 



_B(±L/2) = 0. Similar to the case discussed in section 
|IIIB[ this imposes a strong perturbation on the the low 
energy states as both and \^2) [see Eq. (|4|] have 
to be modified from their bulk form. We study this case 
using the perturbative approach introduced above. The 
unperturbed (71 = 0) states satisfying the boundary con- 
ditions are given by edge states of the form |<&±o) [Eq- 
([18])] (with |Gl), I^l) replaced by [Ofl), 1^^) for right-edge 
states, i.e. X > Q) and the bulk state |$o) [Eq. (12)]. 



Note that the lower energy branch of the edge states, 
elQ{X) = V/2 - uj^y/XiX) (e = R,L), is hole-like and 

crosses the unperturbed bulk level £3*^^ — —V/2. Turn- 
ing on the inter-layer hopping 71 leads to a shift of the 
latter bulk level and its dispersion at the edge, and in 
addition to mixing of the crossing levels and an opening 
of a gap. As in to section [Til B[ we first evaluate the dis- 
persive energy bands £1 _|_(X), £| _|_(X) resulting due to 
mixing with the higher LL n = ±1 to leading order in 
71. The n — ±.1 states of the uncoupled layers are given 
in this case by 



1$ 



±1/ 



1 



± 



V 



1 



) J 



±uj,^Xi{X) , 



(34) 



± 








V 



4i = -^±^cVl + A(X) 



(35) 



where |l'e)|,=±L/2 = 0, |0^) = a|i:,) and A(X), Xi{X) are 
the same as in section |IIIB The resulting perturbative 
expressions for the edge bands dispersing from £1, £2 are 



V 



£ 



2^l\{UQ\^{y±u,.JX{X) 



1,± 



{X) « -±a.,yA(X) ^— ^ (36) 

^ M^M^X{X)[l + X{X)] (a;2 zp 2VuJc^/X{X 



V 2jfvm\0)\' 
2 ^ (A/;')V^A?(X) 
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(37) 



The band elC^) e xhibit s the same behavior as S2{X) 
obtained in section III A [see Eq. ([l6])], which arises in 
both cases from the dominant boundary condition on the 
component B. This corresponds to a moderate hole-like 
dispersion, which interpolates between the bulk energy S2 
and ~V/2 as X is pushed farther and beyond the edge. 
From Eq. (36), the lower branch ef _{X) is also hole- 
like and disperses more steeply. As a result, £|(X) and 
el _{X) tend to cross at X satisfying S2{X) = el _(X). 
As in the case discussed in section |III B[ this crossing 
become avoided and a gap is opening, given (to leading 
order in 71) by 



Ae(^) 



7i 



1(014)1 



(38) 



The resulting edge spectrum is characterized by two sepa- 
rate hole-like bands: one interpolating between the bulk 
state ei — V/2 and a saturated value and one 

starting at £2 and steeply dispersing downwards without 
bound. On top of these, the branch el ^{X) is largely 
particle-like and steeply disperses upward for X near or 
beyond the edge. This behavior is consistent with Fig. 
[Sj It should be noted that the above analysis, based on 
a perturbative expansion in 71, appears to be qualita- 
tively valid even if 71 is not small. As we have argued 



in sections III A and III B the perturbative expansion in 
fact becomes increasingly more justified as X is pushed 
farther over the edge. 



spin-filtering2^, and multicomponent Luttinger liquids. 
These possibilities will be explored in future research. 
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IV. CONCLUSION 



In this paper we have studied edge states of bilayer 
graphene systems in the quantum Hall regime. Our re- 
sults show that a variety of edge state energy structures 
are possible depending on precise boundary conditions. 
In some cases we found that for a continuum model, edge 
states can disperse from a bulk energy value ±V/2 to 
^V/2, while in other cases they may disperse to ±00. In 
yet other cases the edge states may not disperse at all. All 
these behaviors could be understood qualitatively within 
the framework of perturbation theory, and in the first 
of these cases a variational approach allows us to relate 
the edge state dispersion to the problem of edge states 
in single-layer graphene and to the edge dispersion of 
conventional quantum Hall states. The complicated dis- 
persions discussed in this paper yield a variety of possible 
crossings and anticrossings, particularly when spin is in- 
cluded as a degree of freedom and the efi^ects of Zeeman 
coupling are considered. This rich set of possible spectra 
for the edge states of bilayer graphene in a magnetic field 
suggest a variety of possibilities for physical phenomena 
at the edge, including counterpropagating edge states. 
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